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Abstract: We investigate the weak electric and magnetic dipole moments of top quark- 
Z boson interactions at the Large Hadron Collider (LHC) and the International Linear 
Collider (ILC). Their vanishingly small magnitude in the Standard Model makes these 
couplings ideal for probing New Physics interactions and for exploring the role of top quarks 
in electroweak symmetry breaking. In our analysis, we consider the production of two top 
quarks in association with a Z boson at the LHC, and top quark pairs mediated by neutral 
gauge bosons at the ILC. These processes yield direct sensitivity to top quark-Z boson 
interactions and complement indirect constraints from electroweak precision data. Our 
computation is accurate to next-to-leading order in QCD, we include the full decay chain 
of top quarks and the Z boson, and account for theoretical uncertainties in our constraints. 
We find that LHC experiments will soon be able to probe weak dipole moments for the 
first time. 
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1. Introduction 

Twenty years after the first observation of the top quark, experiments at the Tevatron 

and the LHC have provided us with a detailed picture of top quark properties and inter¬ 
actions. Today, we have precise knowledge of its basic properties such as the mass, the 
pair production cross section, spin correlations and the kT-helicity fractions, all of which 
are highly consistent with the Standard Model (SM) predictions. Moreover, the produc¬ 
tion of tt pairs at the Tevatron and the 7 and 8 TeV LHC spans a wide range of energies 
which profoundly probes the QCD production mechanism in different partonic channels, 
as well as the charged-current decay dynamics. This has resulted not only in a detailed 
understanding of the top quark, but has also led to strong exclusion limits for New Physics 
scenarios. 

One aspect of top quark physics which is far less explored are interactions with neutral 
electroweak gauge bosons and the Higgs boson. Sensitivity to these couplings arises in 
hadronic collisions through the partonic subprocess qq —> Z* /'y* —> tt and through higher- 
order electroweak corrections to pp —> tt. However, at the Tevatron and the LHC the 

Z*/ 7 * mediated process is overwhelmed by the strong production mechanism, and event 
numbers at high energies are too low to facilitate a coupling analysis via electroweak loop 
corrections. A more promising avenue for anomalous coupling studies is the investigation of 
the processes pp —tt + V with V = Z, 7 , H which yield direct sensitivity without intrinsic 
dilution by QCD effects or loop suppression. Their relatively high production threshold of 
2 nit + W requires a sufficiently large center-of-mass energy which is, however, adequately 
reached by the LHC. ATLAS and CMS have already measured the tty cross section at 
Run I [1,2], and there is evidence for the production of the ttZ final state [3,4], Run 
II promises 10 1 * * * 5 ttZ events per 100 fb _1 luminosity, and about half as many ttH events. 
Hence, the study of tt + V final states at the LHC will soon open a new chapter in top 
quark physics and allow anomalous top quark couplings to electroweak gauge bosons to be 
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probed directly for the first time. Ultimately these efforts will provide new insights into 
the role of top quarks in electroweak symmetry breaking and the effects of New Physics 
dynamics. 

Current and past collider experiments already provide some information on anomalous 
top quark couplings. For example, electroweak precision data from LEP experiments can be 
cast into constraints on the parameters e\, £ 2 , £3 and £b [5-7], all of which are sensitive to 
top quark anomalous couplings. The authors of Ref. [ 8 ] use these parameters to constrain 
the weak magnetic dipole coupling of the top quark. Flavor changing penguin transitions 
measured at R-factories and kaon experiments are also a sensitive probe of top quark in¬ 
teractions. The authors of Ref. [9] use rare B s and K meson decays to constrain vector 
and axial ttZ couplings, whereas in Refs. [10,11] B decays and asymmetry observables are 
used to put limits on the electromagnetic dipole moments. To the best of our knowledge, 
there are at present no constraints on the weak electric dipole moment. Although con¬ 
straints from LEP and heavy meson experiments can be quite strong, it should be noted 
that they often arise in a convoluted and rather indirect way. For example, the analysis 
of Ref. [ 8 ] shows that constraints involving the top quark weak magnetic dipole moment 
involve three other anomalous parameters, leading to an ambiguous interpretation. This is 
not surprising given that top quarks were never observed at LEP and sensitivity only arises 
through virtual top quarks in electroweak loop corrections. Similarly, in Z boson medi¬ 
ated penguin transitions the top quarks only enter through virtual exchange and involve 
additional couplings which can be anomalous. In contrast, sensitivity in ttZ production at 
the LHC arises directly through interactions of on-shell particles without much dilution by 
other possible anomalous couplings. 

Looking further into the future, the prospect of an International Lepton Collider (ILC) 
offers the possibility of extending the LHC top quark program [12]. The clean collider set¬ 
ting in e + e _ collisions and the option for polarized beams presents an optimal environment 
for precision studies of SM processes. In particular, if the center-of-mass energy is above the 
top-pair production threshold, the couplings of the top quark to the neutral gauge bosons 
are directly accessible through the process e + e _ Z* 1 7 * —> tt. At energies far above the 
tt threshold even the process e + e~ —>• tt + H opens up and allows precision studies of the 
top quark Yukawa coupling. Moreover, the study of anomalous couplings in e + e _ collisions 
goes beyond the exploration of New Physics models, as it is currently also important for 
estimating the capabilities of a future linear collider. Reliable and accurate benchmarks 
for how well an ILC experiment can pin down fundamental interactions is a major criterion 
for future planning that needs to be addressed decades before construction and data taking. 

In this paper we focus on a particular kind of anomalous couplings: the weak mag¬ 
netic and electric dipole moments of top quark-Z boson interactions. In the SM those 
dipole moments are zero at tree level and higher-order electroweak corrections introduce 
values which are too small to be measured anytime soon. Hence the study of anomalous 
dipole moments provides an ideal place to look for New Physics in the third generation 
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electroweak sector. Various Beyond the Standard Model scenarios which address the mech¬ 
anism of electroweak symmetry breaking predict sizable dipole moment couplings, as well 
as deviations of the vector and axial ttZ couplings from their SM values [13-19]. This 
raises the question of how well collider experiments can access their parameter space. This 
problem has been investigated in Refs. [20,21] for the LHC, and in Refs. [22,23] for the ILC. 
The former analysis clearly indicates that the leading-order (LO) scale uncertainties are the 
biggest limiting factor in setting more stringent constraints at the LHC. We recently began 
to address this problem by constraining the vector and axial ttZ couplings at the LHC 
through next-to-leading order (NLO) in perturbative QCD [24], finding an improvement 
by as much as a factor of two compared to a LO analysis. In this paper we expand our 
efforts to include the anomalous weak dipole interactions at the LHC, and also consider 
constraints on the ttZ couplings at a future ILC with \/s = 500 GeV, both at NLO QCD 
accuracy. We also briefly comment on limits achievable in 100 TeV pp collisions. 

The remainder of the paper is organized as follows. In section 2, we present precise 
definitions of the ttZ couplings, briefly discussing their relationship to operators in an 
effective Held theory , and we outline the details of the calculation. We then consider LHC 
constraints in section 3 and ILC constraints in section 4. We conclude in section 5. 


2. Anomalous top quark couplings 


We begin with the Standard Model Lagrangian for ttZ interactions 


C^ = eu{p t ) 


u( n z, SM 

r( c v +75 l a ) 


v(pt)Z, 


/M 


( 2 . 1 ) 


where e is the electromagnetic coupling constant. The vector and axial couplings are 
T 3 — 2 Qt sin 2 


a 


Z,SM 


V 


2 sin 9 m cos 9,, 


0.24, 


C 


Z,SM 


~T t 3 

2 sin 9 W cos 9,, 


-0.60, 


( 2 . 2 ) 


with Qt = 2/3, T 3 = 1/2, and the weak mixing angle 9 W . To parameterize possible effects 
from physics beyond the SM the Lagrangian can be extended to 


Aiz = eu {p t ) 


Ppiv + 75 Cy A ) 7 ip 2 ,v + 1756 ^ 4 ) 


M z 


v (Pt)Z) 


/M 


(2.3) 


with = 5 [ 7 ^, 7 ^] and q v = (pt —Pt)v- The coupling C% v corresponds to the weak mag¬ 
netic dipole moment and Cp 4 to the CP-violating weak electric dipole moment. Higher- 
order corrections within the SM effectively induce finite values for these dipole moments. 
However, their size is very small, in particular C%y ~ 10 -4 [25] and only receives 
contributions at three-loops [13,26]. The anomalous couplings are conveniently parame¬ 
terized by operators in an effective field theory (EFT) that respects the symmetries of the 
SM. Following Ref. [27], the four couplings C/y 2 V , A in Eq. (2.3) have a simple translation 
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into four dimension-six operators 


Cf >v = 0 M + (£) Re [2Cg 33) - 01] , Cf A = Cf* + 


A 2 


Re 


oai(3,33) . ^33 


(2.4) 


C%y — V2 


A 7 


Re 


C W0W - Sw0Bcj> ■ OfA — ( -^2 ) Irn [cyr'C^y - sw/C 13 ^] 


Hence, any constraint on the anomalous couplings can be translated into constraints on a 
combination of effective operators. After UV renormalization and evolution of the operators 
from the EFT scale A down to the electroweak scale, the operators mix at NLO and 
additional logarithms log (Mg/A) appear. In this paper we refrain from a comprehensive 
operator analysis and present constraints on the anomalous couplings only. A forthcoming 
publication on anomalous tty couplings will allow us to come back to this issue in a more 
general way. 

In Ref. [24] we studied the anomalous couplings C0 , A in the process pp —> ttZ and 
treated the Z boson in the narrow-width approximation (although spin-correlated decays 
into leptons were fully included). In this work, we will investigate the V / A couplings 
which involve a momentum-dependent term <j ,w q v Z tI , which is to some extent sensitive to 
off-shell effects of the Z boson. We therefore extend our previous calculation and allow 
the Z boson to be off-shell, and include intermediate off-shell photons decaying into lep¬ 
tons. Consequently, picking up interference between Z and 7 contributions implies that a 
possible deviation from the SM couplings can no longer be unambiguously traced back to 
anomalous ttZ interactions. This issue is also apparent in the effective field theory where 
anomalous ttZ and tty interactions share certain effective operators (cfg. Eq. (39) and Eq. 
(41) in Ref. [27]) and can lead to correlated shifts. We will show in the next section that 
the contribution from intermediate photons in ttZ final states at the LHC is negligible as 
long as the invariant mass of the two leptons is not extremely far away from Mz- Hence, 
for the scope of this work, we use the SM tty coupling, and treat ttZ production at the 
LHC as a pure and direct probe of top-Z interactions, even though intermediate photons 
are present. At the ILC, the situation is different. The intermediate off-shell Z boson and 
photon mix at the fixed collider energy and thus the contribution from the intermediate 
photons is not negligible. Possible avenues to resolving the ambiguity have been explored in 
Refs. [22,23], facilitating polarized e + e _ beams. However, in our study of ILC constraints 
in section 4 we will not address this issue, and will only consider anomalous ttZ couplings, 
again keeping the tty couplings at their SM value. We will, however, come back to this 
issue in a separate publication on anomalous pp —> tty production at the LHC. 


On the technical side, our calculation builds upon the results of Ref. [24], We calculated 
the NLO QCD corrections to pp —>• ttZ with leptonic branchings of the Z boson and include 
top quarks decays with NLO spin correlations. The top quarks are treated in the narrow- 
width approximation which is parametric in Tt/mt as long as the invariant mass of the 
top quark is integrated over. Similar calculations for on-shell tops and parton showered 
decays have been presented in Refs. [28,29], For this study, we extend our framework by 
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implementing the weak dipole interactions in Eq. (2.3), and by adding a full simulation of 
the process e + e - —> tt —> bjj M~i> at NLO QCD, allowing for all possible ttZ anomalous 
couplings. In this context we briefly mention that the weak dipole couplings in Eq. (2.3) 
introduce additional UV-singularities at NLO QCD. We subtract these divergences through 
a redefinition of the C% v / A couplings 


u 2,V/A 


r z 

u 2 y/A 


1 + 


Mren 

M| 


E(l + e) 


(2.5) 


and confirm their universality in the processes pp -+ ttZ and e + e - tt. In the EFT this 
procedure is equivalent to an operator renormalization which arises through the mixing of 
operators at 0(a s ). 


3. LHC Constraints 

In this section, we consider the process pp —> tt Z /y* —> bjj M~ v £ - £ + with one W boson 
decaying hadronically, the other decaying leptonically, and a leptonic branching of the Z 
boson. Thus, the final state consists of three charged leptons, four jets (two of which 
may be fr-tagged), and missing energy. In our numerical results we sum over the possible 
combinations of the charged leptons = e^, . We concern ourselves exclusively with 

the LHC operating at a center-of-mass energy y/s = 13 TeV and fix the input parameters 
as follows: 


mt = 173 GeV, 

M z = 91.1876 GeV, 

C F = 1.166379 x 10" 5 GeV -2 , 
T^° = 1.4957 GeV, 

T^ = 2.0455 GeV, 


mb = 0 GeV, 

M w = 80.385 GeV, 

T z = 2.4952 GeV, (3.1) 

r t NLO = 1.3693 GeV, 
r^ LO = 2.1145 GeV. 


We use the central factorization and renormalization scale po = mt+Mz /2 and estimate the 
scale uncertainty by varying it by a factor of two in either direction. Parton distributions 
are taken from MSTW2008 [30] with a s (Mz ) = 0.13939 and a s (Mz ) = 0.12018 at LO 
and NLO QCD, which we evolve to the renormalization scale using one-loop and two-loop 
running, respectively. We choose acceptance cuts on the final state to correspond closely 
to those used in a CMS analysis of the full 8 TeV dataset [3]: 


> 20 GeV, 
p j T > 30 GeV, 
Rij > 0.3, 

mu — Mz\ < 20 GeV. 


\y e \ < 2.4, 
W\ < 2-4, 
Ru > 0.3, 


The jets are defined according to the anti-fcx algorithm [31] with R = 0.5. With this setup, 
we obtain a LO and NLO QCD cross-section of 


LO _ 2 Q4+36% 
a - z.uu_ 24% 


fb, 


NLO _ q qo _ l _ 19% 

a - 3.38_ lg% 


fb, 


(3.3) 
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Figure 1: NLO QCD distributions of A <j>u and Pt,z for the 13 TeV LHC, using the cuts of 
Eq. (3.2). Shown are the results allowing the Z boson to be off-shell (a), keeping the Z boson 
on-shell (b), and allowing the Z boson to be off-shell but removing the off-shell photon coupling to 
the top quarks (c). The lower pane shows the ratios (b)/(a) and (c)/(a). 


where the lepton multiplicity of eight has been included. The scale uncertainty is obtained 
by varying the scale by a factor of two in either direction. The iC-factor of 1.66 is slightly 
larger than the one found in Ref. [24] for the tri-leptonic final state. This is explained by 
the more stringent jet cuts used in the present analysis, which are passed more readily 
when an additional jet is present in NLO kinematics. 

We begin the discussion of results by examining the impact of the Z boson off-shellness 
and the photon contribution and compare it to the situation in which the Z boson is treated 
in the narrow-width approximation. In figure 1, we show two distributions: the azimuthal 
angle between the two leptons originating from the Z boson decay, A (j>u, and the transverse 
momentum of the Z boson, pt,z- Allowing the Z boson to go off-shell decreases the cross 
section by about 3% at NLO with very little shape change in either distribution. This is 
fully consistent with the parametric estimate Tz/'Mz ~ 0.03 as obtained from the narrow- 
width approximation. Removing the contribution of the photon results in a further uniform 
decrease by 1%. This small effect is easily explained by the window cut on the invariant 
mass of the lepton system, given in Eq. (3.2), which forces the photon to be far off-shell. 
This implies that our setup has very little sensitivity to the couplings of the top quark to 
the photon, and we therefore keep tt 'y couplings at their SM values without introducing 
anomalous interactions. 

We now consider the impact of anomalous weak dipole couplings C z v and (7^ A on the 
ttZ results. For illustration, we choose C z v = C Z A = 0.2, but keep C z v and Cf A at their 
SM values, and find 

£x LO (Cfv = C 2 Z A = 0.2) = 3.51 fb, <t nlo (C^ v = C 2 Z A = 0.2) = 5.99 fb. (3.4) 

These cross sections are notably larger than the corresponding SM ones in Eq. (3.3). More¬ 
over, this choice of anomalous couplings leads to substantial shape changes, as shown in 
figure 2. We present the normalized A^ and pt,z distributions and compare anomalous 
versus SM coupling hypotheses. The shape changes are significantly stronger than those 
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Figure 2: NLO distributions of A<ji« and pt,z for SM tiZ couplings and with anomalous dipole 
couplings C. f v = Ci A = 0.2. The distributions are normalized to the overall cross section. The 
cuts of Eq. (3.2) are applied. 


obtained from anomalous Cf v and Cf A couplings as presented in figure 5b of Ref. [24], 
The reason is clear: varying the values of Cf v or Cf A only changes the relative strength 
of SM vector and axial couplings. In contrast, finite contributions of the dipole couplings 
C.fy or C^ A introduce a new Lorentz structure to the ttZ vertex which is proportional to 
the Z boson momentum vector. This leads to significantly different spin correlations as 
well as enhancements in the tail of energy-related distributions. In the following analysis 
we will use the two observables presented in figure 2 as they provide strong sensitivity to 
the anomalous couplings 1 . The steady enhancement of anomalous cross sections in the 
Z boson transverse momentum distribution will violate unitarity at very high energies. 
This is related to the transition when the dipole couplings turn into g 2 -dependent form 
factors, or equivalently, when the EFT breaks down at energies ~ Aeft- One solution 
is to introduce a damping factor (1 + q 2 / A 2 ) -2 [20] which prevents unitarity violation at 
high energies. In our analysis we simply choose to disregard all normalization and shape 
information at pt,z > 500 GeV. Assuming that BSM physics enters at scales far above 
500 GeV this should be a safe procedure. 

We now proceed to studying the constraints that LHC experiments will be able to 
place on the ttZ couplings with 13 TeV data from Run II. For reasons of completeness 
and comparability with the ILC results in section 4, we first repeat our earlier analysis [24] 
by independently varying Cf v and A while keeping y and Of A zero. In fact, we 
marginally extend the framework of Ref. [24] by accounting for off-shell Z bosons and 
photons which were neglected earlier. We then vary C^y and C^ A independently while 
keeping Cf v and Cf A at their SM values. Thus we will obtain constraints in the parameter 
spaces of (Cf v , Cf A ) and {C^y, A )- Our strategy closely follows that of Ref. [24]: We 

construct a fitting function to generate LO and NLO kinematic distributions for a large 

1 We note that A <j>u and pr,z are by no means the only observables whose shape has strong discriminating 

power. For example, the angle between the semi-leptonically decaying top and the charged lepton originating 
from this decay in the top quark rest frame is another ideal observable. 
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Figure 3: Significance of deviations from the SM for the couplings C-f v and Cf A with 30, 300, 
and 3000 fb^ 1 of data at the 13 TeV LHC. The results are obtained from rate and A (f>u shape 
information using LO predictions shown on the left and NLO QCD predictions shown on the right. 


number of anomalous couplings. We then use a log-likelihood ratio test to determine the 
significance with which non-SM values of the couplings can be separated from the SM hy¬ 
pothesis, assuming that the SM hypothesis is true. In the event of the experimental results 
being compatible with the SM, this can be interpreted as the approximate significance 
with which a given value of the ttZ couplings can be excluded. We include the effects of a 
theoretical uncertainty Atheor. unc. hr the likelihood analysis through the overall rescaling of 
the results under both hypotheses, such that the difference between the two cross sections 
is minimized. The value of the theoretical uncertainty is taken to be A t heor. unc. = 30% at 
LO and A t h eor . unc . = 15% at NLO, corresponding roughly to the scale uncertainties from 
missing higher order corrections. 
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The constraints on C± v and Cf A . obtained using the A (f)^ distribution, are shown in 
figure 3 for an integrated luminosity of 30, 300, and 3000 fb” 1 . The area outside of the 
blue area roughly corresponds to couplings that can be excluded at 68% confidence level 
(C.L.); the area outside the red area can be excluded at 95% C.L. Thus we expect coupling 
constraints, using 300 fb -1 of data and employing NLO QCD results, of approximately 
-0.75 < Cf A = -0.60 < -0.5 and -0.7 < Cf y = 0.24 < 0.7 at the 95% C.L. The 
improvement when using NLO QCD over LO predictions is immediately apparent from 
the figures. This is mainly a consequence of the reduced scale uncertainty and the larger 
event rate due to a positive perturbative correction. 

Turning now to the dipole moment couplings C^y and Ctf A , we focus on the constraints 
obtained from the pt,z distribution in figure 4. We choose a different color scheme to 
that of figure 3 to emphasize the different anomalous couplings. We note that our LO 
constraints are consistent with those of Baur et. al. [20], and including NLO corrections 
improves the constraints by 20%-40%. At this accuracy, the LHC experiments will be 
able to probe values of |C’^ V/A <0.3 with 30 fb” 1 of data, improving to |C^ V , A | <0.15 
with the full 300 fb -1 dataset and reaching | ~ 0-08 with the high-luminosity 

upgrade to 3000 fb” 1 of data. The effect of the NLO corrections is slightly smaller than 
that observed for the constraints on V / A - This can be understood from our earlier 
observation (figure 2) that the shapes of the distributions have greater sensitivity to dipole 
moment couplings. Although the overall rate is sensitive to variations in the couplings, 
the shape information dominates and provides the main contribution to our constraints. 
Our treatment of the theoretical uncertainty involves rescaling the overall rate, so that 
reducing the scale uncertainty has the effect of substantially improving the contribution to 
the constraints coming from the rate alone, while having less of an effect on the constraints 
that are obtained from the shape. Since our constraints are dominated by the latter, the 
improvement in the constraints when the scale uncertainty is decreased is mild. This also 
explains why we see a greater improvement at lower luminosities, as this is where the rate 
alone is more important in obtaining the constraints. 


Aside from a high-luminosity LHC and an ILC, the particle physics community is 
also considering the construction of a Future Circular Collider (FCC) with proton-proton 
collisions at a center-of-mass energy of 100 TeV [32-34]. While this machine is usually 
thought of in the context of exploring a new energy regime, SM particles will be produced 
in great abundance and it is interesting to know how analyses compare to those from the 
clean ILC environment. For our purposes, the pp ttZ cross section at the FCC is about 
50 times larger than at the 13 TeV LHC, resulting in over 10 8 ttZ events expected with 
10 ab” 1 of data. In a rough estimate we repeat our analysis for the FCC at leading-order, 
assuming that the theoretical uncertainty can be reduced to A t heor. unc. = 5%. Our findings 
indicate that FCC constraints improve by factors of 3 — 10 compared to the 3 ab” 1 LHC. 
In relation to our ILC results (presented in the following section) the FCC constraints are 
stronger by factors of a few. 
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Figure 4: Significance of deviations from the SM for the dipole couplings Cf v and Cf A with 30, 
300, and 3000 fb -1 of data at the 13 TeV LHC. The results are obtained from rate and pt,z shape 
information using LO predictions shown on the left and NLO QCD predictions shown on the right. 


4. ILC Constraints 

We now turn our attention to a future e + e _ collider. Such a machine offers several ad¬ 
vantages over a hadron collider, especially in performing SM precision measurements [12]. 
It provides an experimentally clean environment without hadronic activity in the initial 
state, the collision energy is accurately known, and the possibility of polarized beams of¬ 
fers additional handles on analyses. At present, the most advanced proposal for an e + e - 
collider is the International Linear Collider (ILC) [35] with a design center-of-mass energy 
of \/s = 500 GeV and the possibility for an upgrade to 1 TeV. Precision top quark studies 
form an integral part of the physics program for this machine since at y / s = 500 GeV top 


10 





































pairs are produced numerously well above threshold. Moreover, the production proceeds 
through an intermediate Z/'y exchange which makes this process sensitive to electroweak 
top quark couplings already at leading order. 

The ability of an ILC experiment to measure these couplings was analyzed in Refs. [22, 
23], using as discriminating observables the cross section, the scattering angle, the forward- 
backward asymmetry, and the slope of the helicity angle asymmetry. These works also take 
advantage of two beam polarizations to disentangle the tt'y and ttZ couplings, and find 
sensitivities for the top quark couplings between 2-6 permille (depending on the coupling) 
assuming 500 fb” 1 of data. In this paper we use a slightly different approach to estimate 
our constraints. Adopting the same collider parameters of \fs = 500 GeV and 500 fb -1 
of data, we assume a simplified experimental procedure with unpolarized beams and only 
one single discriminating observable. On the other hand, we account for higher-order 
QCD corrections and include the corresponding theoretical uncertainties in our constraints. 
Hence, this procedure matches the one applied in our LHC analysis in section 3, enabling 
a direct comparison. Following Ref. [22] we use the scattering angle of the top quark to 
probe anomalous couplings. We verified that this observable has a strong sensitivity to the 
ttZ couplings and we consider the lepton+jet final state which allows a full reconstruction 
of the top quark momenta. We apply the following cuts to the tt final state 

p ( T > 20 GeV, |/| < 2.5, 

p{ > 20 GeV, \y j \ < 2.5, (4.1) 

p™ ss > 20 GeV, RjikJ 1 ) = 0.4 

and find cr NLO = 90.0 ±0.8 fb. This promises 45000 top quark pairs from 500 fb -1 of data. 
Our prediction includes NLO QCD corrections in production and decay of the top quarks 
maintaining full spin correlations within the narrow-width approximation. As pointed out 
in section 2, the anomalous dipole couplings in e + e” —> tt are UV-renormalized in exactly 
the same way as in hadronic collisions, given by Eq. (2.5). The QCD NLO corrections to 
top-pair production in e + e” annihilation have been known for a long time [36-40]. NNLO 
QCD corrections to the top-pair production at threshold have received significant attention 
in the past (see [41] and reference therein), and progress is being made on the N 3 LO 
corrections [42-44] at threshold. In the continuum, certain differential results at NNLO 
QCD have recently been presented recently [45-47]. Both the NLO scale uncertainty and 
the NNLO corrections at ^/s = 500 GeV are at the level of l%-2%, indicating a very good 
perturbative convergence [45]. NLO electroweak corrections to the overall cross section are 
at the level of a few percent, but Refs. [48,49] report shape changes of up to 10% in the 
cos 6 differential distribution. 

An important part of our analysis is the inclusion of theoretical uncertainties in the 
constraints. Since we neglect known electroweak corrections we will present two scenarios: 
(1) We assume a theoretical uncertainty of as much as 10% in our log-likelihood ratio test to 
account for the missing electroweak corrections in our predictions; (2) we estimate future 
results with 2% uncertainty assuming that NNLO QCD and 0(a ) corrections including 
anomalous couplings are known. Figure 5 presents the results for the first scenario. Sensi- 
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Figure 5: Sensitivities of deviations from the SM for the couplings Cf v and C£ A (left), and C£ v 
and C^ A (right) for the 500 GeV ILC with 500 fb -1 of data. The constraints assume NLO QCD 
predictions and a total uncertainty of 10%. For ease of comparison the scale of the large panes is 
kept the same as in our LHC analysis. 


tivity to deviations of Cf v and (%f A is shown on the left, and the sensitivity to deviations 
°f and C,f A is shown on the right. For ease of comparison the scale of the plots 
is kept the same as for the LHC (figure 3,4). Exactly the same is shown in figure 6 for 
scenario (2) with 2% theoretical uncertainty. In view of future results we choose to discuss 
the latter. We find that vector and axial couplings can be constrained to approximately 
Cfv = 0.241 q°® and Cf A = — 0.60^o 04 the 95% C.L. Comparing these with the con¬ 
straints presented in section 3 for the high luminosity LHC, we see an improvement by a 
factor of three on the Cf A constraints, and an improvement by a factor of approximately 
six on Cfy At both the LHC and the ILC, we see a greater sensitivity to Cf A than to 
Cfy We also note that fixing one of the two couplings to its SM value improves the ILC 
constraints on the other coupling by a factor of about three. 

Turning now to the weak dipole moments (figure 6, right), we find the constraints 
|Cf A I 1$ 0.08, and —0.02 < Cf v 1$ 0.04 at the 95% C.L. Comparing to the LHC constraints 
presented in section 3, we see that the constraints on Cf A from the ILC are comparable 
to those from the LHC with 3000 fb -1 of data, as is the upper limit on Cfy However, 
the asymmetric pattern of the ILC constraints on Cf v allows an improvement on the 
lower limit of this coupling by a factor of three compared to the LHC. We emphasize 
that the constraints presented here are purely a statistical comparison between different 
calculations, and do not take into account experimental effects, which in general will weaken 
them. Since we expect the experimental effects to be larger at the LHC than the ILC, the 
above discussion should be viewed as a conservative estimate of the relative improvement 
of the ILC constraints over the LHC constraints. 

It is worthwhile to compare the results obtained here with those of Amjad et. al. [23], 
who also consider the ILC running at yjs = 500 GeV and with an integrated luminosity of 
500 fb -1 . Varying one coupling at a time, they find constraints at the 68% C.L. of 4 per 
mille for Cfy. 6 per mille for C, and about 1 per mille for Clfy They do not constrain 
the C*^ A coupling. We also note that, although Ref. [23] contains a thorough exploration 
of the experimental aspects, the quoted sensitivities do not include systematic errors. If we 
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Figure 6: As for figure 5, but assuming a theoretical uncertainty of 2%. For ease of comparison 
the scale of the large panes is kept the same as in our LHC analysis. 


assume SM values for all but one coupling, and use A t h e or. unc. = 2%, then at the same C.L. 
we find similar constraints on Cf A . However, our constraints are a factor of five weaker 
for Cf v and Cf v . This can be understood by several differences in the analyses: In this 
paper we only use one single differential observable and assume unpolarized beams, and, 
furthermore, we include theoretical uncertainties from scale dependence in our constraints. 
We stress that this last factor is an unavoidable source of systematic uncertainty which 
should be taken into account when projecting future constraints. Nevertheless it is clear 
that a strong improvement can be obtained once polarized beams and a multi-variate 
analysis are employed. 

5. Conclusions 

In this paper we have continued our work on directly constraining ttZ couplings in collider 
experiments. We discuss the study of weak magnetic and electric dipole moments, two 
couplings which are an excellent probe of new physics given their small SM values. In 
our analysis we consider anomalous production of top quark pairs in association with a Z 
boson at NLO QCD and generate various differential distributions over large grids in the 
anomalous coupling space. We use the obtained shape and normalization information in a 
log-likelihood ratio test which includes leading theoretical uncertainties in order to reliably 
estimate future constraints. 

We find that experiments at the 13 TeV LHC will soon be able to explore weak dipole 
couplings in ttZ production, and that including higher-order predictions boosts the range 
of exclusion limits. For example, the magnetic and electric tiZ dipole moments can be 
constrained to ±0.15 with 300 fb -1 in the tri-lepton final state. This is an improvement by 
25% compared to an analysis at leading order. A further improvement to ±0.08 is possible 
once 3000 fb” 1 of data is accumulated. Similar constraints can be extracted from LEP 
electroweak precision data, and possibly from future H-physics analyses. However these 
indirect limits are unavoidably diluted by assumptions on other anomalous contributions 
which enter the observables. On-shell production of ttZ at the LHC, on the other hand, is 
a much more direct probe and will soon yield unambiguous and stringent results. 
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We also estimate constraints from an NLO QCD calculation of ti production at the 
ILC, and find constraints of —0.64 < Cf A < —0.56 and 0.13 < < 0.32, approximately 

a factor of three to five better than those obtainable by the LHC. For the weak dipole 
couplings we find —0.08 < C^ A < 0.08 and —0.02 < Clfy < 0.04, which represents a factor 
of three improvement for C^y 

Finally, we note that if anomalous ttZ couplings are expressed through dimension- 
six operators of an effective theory they share operators with tt'y interactions. Since the 
process pp — > tt'y is accessible at the LHC this offers the unique opportunity to further 
constrain these operators and pin down top quark electroweak interactions. Similarly, ILC 
experiments will be able to disentangle Z and 7 interactions to top quarks if polarized 
beams are available. 
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